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1. Introduction 

Multi-matrix integrals of various types appear in many mathematical and physical 
applications, such as combinatorics of graphs, topology, integrable systems, string theory, 
theory of mesoscopic systems or statistical mechanics on random surfaces. 

A general Q-matrix integral of the form 



usually goes over the N x N hermitian, real symmetric or symplectic matrices M q with 
the action S and the measure symmetric under the simultaneous group rotation: M q — > 
Q + M q Q. Some other multi-matrix integrals, such as these with complex matrices or with 
general real matrices, can be reduces to those three basic cases. 

We will consider here only the case of hermitean matrices for which O belongs to the 
U (N) -group. 

In many applications "to solve" the corresponding matrix model usually means to 
reduce the number of variables by explicit integrations over most of the variables in such 
a way that instead of QN 2 original integrations (matrix elements) one would be left in 
the large iV limit only with ~ iV integration variables. In this case the integration over 
the rest of the variables can be performed, at least in the widely used large N limit, by 
means of the saddle point approximation. A more sophisticated double scaling limit [|1| is 
also possible (if possible at all) only after such a reduction. The key of success is in the 
fact that after reduction the effective action at the saddle point is still of the order ~ N 2 
whereas the corrections given by the logarithm of determinant of the second variation of 
the action cannot be bigger than ~ N (the "entropy" of the remaining variables). The 
problem is thus reduced to the solution of the "classical" saddle point equations, instead 
of the "quantum" problem of functional (in the large N limit) integration over the original 
matrix variables. 




5=1 



(1.1) 
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Such an explicit reduction of the number of " degrees of freedom" is in general possible 
only for a few rather restricted, though physically and mathematically interesting, classes 
of multi-matrix integrals. The purpose of our present notes is to review the basic old and 
new methods of such a reduction. Before going to the particular cases let us stress the 
importance of the search for new methods of such a reduction: any nontrivial finding on 
this way leads immediately to numerous fruitful applications. 

2. Some old examples 

The best known example of such a reduction of the number of degrees of freedom is 
the one matrix integral: 

Z = J d N2 MexpNTrS(M) (2.1) 
where S(M) is an arbitrary function of one variable. Let us use the decomposition: 

M = O+xO (2.2) 

where x = diag(x\, • ■ • , xn) is a diagonal matrix of the eigenvalues and O is the U(N) 
group variable. The corresponding (Dyson) measure can be written as: 

N 

d N M = d[fl} u{N) A 2 (x) Y[ dx k (2.3) 

fe=i 

where A(x) = Yli>j( x i ~ x j) * s the Van-der-Monde determinant. The integrand as an 
invariant function does not depend at all on O (the integration over it produces just a group 
volume factor which we will always omit). The remaining integral over the eigenvalues 
reads: 

r N 

Z= l[dx k exp[NS(x k )}A 2 (x) (2.4) 
J k=i 

In the large N limit the corresponding saddle point equation takes the form 
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= S'(x k ) + ±-Y = (2.5) 



These arguments were successfully used for an interesting combinatorial problem: enumer- 
ation of graphs of fixed two dimensional topologies M, B. There exist powerful methods 
to analyze this equation but it is not our present goal to review them here. 
The next fruitful example is the so called two matrix model: 

Z = J d N2 A d N2 B exp iVTr ( - A 2 - B 2 + cAB + U{A) + V(B)) (2.6) 

where U and V are some arbitrary functions of one variable. After the decomposition 
A = Q^x^2: B = 0^"y02 we are left, due to the term Tr(AB) in the action, with one 
nontrivial unitary integral over the variable O = fiifij. Fortunately, this integral was 
explicitly calculated by Charish- Chandra, Itzykson and Zuber 0,0: 

d[n] u(N) e W Tr(n+xny) = J] fc! A( ; J )A(y) (2-7) 



Substituting (|2.7|) and the Dyson measure ( |2.3| ) into ( |2.6| ) we are left again with only 2N 
variables Xk and and we can write again the saddle point equations in the large N limit. 
They are more complicated than in the one matrix integral but can be nevertheless solved 
quite explicitly. The first solution of that kind was found in |5| in an indirect way, using 
the method of orthogonal polynomials, but the direct solution is also possible (see 0). 

This model was used in to solve exactly the first example of new statistical me- 
chanical models of interacting spins on random planar graphs: in this case it was a model 
of Ising spins on random planar graphs. 

An obvious generalization of the two matrix model is the matrix chain model: 

» Q Q Q-i 

Z = jl[d N2 M q e W Tr(J2V q (M q ) + J2 M p-i M p) ( 2 - 8 ) 

One easily notices that the same unitary decomposition M q = Vt^x q Vt q leads to Q — 1 
independent integrals over the variables U q = fi+^Og of the type (f2.7|). We are left 
again with only QN eigenvalues instead of QiV 2 matrix elements and are ready to apply 
the saddle point approximation to this integral. This model was first analyzed by the 



method of orthogonal polynomials by ||. It was shown in || that by special choices of 
the potential V the model can be described by the KP integrable flow with respect to the 
coupling constant of the potential. 

Note that if we imposed the periodicity condition M\ = M q on this matrix chain and 
add the term M\M q to the action the problem would become much more complicated (and 
actually not solved so far), since this would give an extra condition \\ U q = I making the 
variables U q not independent. 

Another solvable matrix chain describing the statistical RSOS RSOS models on ran- 
dom planar graphs was proposed and solved in |TTJ. Similar models were considered in 

m 



Some multi-matrix models can be reduced to the solvable ones by means of simple 
matrix integral transformations. The first example of such transformation was described 
in the paper |12] for the matrix integral describing the Q-state Potts model on random 



dynamical planar graphs. Its partition function is 

Q Q Q 



Z = J f[d N M q expTr(J2v q (M q ) + J2 M p M «) ( 2 ' 9 ) 

q = l q=l P,Q=1 



5=1 q = l p,q= 

One can represent the last factor under the integral as 

f 1 Q 

/ d N X exp Tr ( - -X 2 +X^M q ). 

J q = l 

Let us consider the case V\ = ■ • • = Vq = V. Then the whole integral can be expressed as 
Z = J d N X exp(~^TrX 2 ) j d N *M exp Tr(XM + V(M)) Q (2.10) 



The integrals in (|2.10|) can be reduced to the eigenvalues: in the integral under the power 



the only nontrivial "angular" integration over the relative U(N)- n angle" can be done by 
means of the formula (|2.7|) and the external one will also depend only on the eigenvalues 



of X. The solution of the corresponding saddle point equations was found in [13] and 



analyzed in and [ |15|| . 

Combining these methods in the obvious ways one can generalize the large iV solv- 
ability on a certain larger class of multi-matrix models. 



3. Matrix Quantum Mechanics 

In the limit when Q — > oo and with the special scaling of coupling constants the 
matrix chain (|2.8| ) becomes matrix quantum mechanics. It is defined by the Hamiltonian 

H M = -A M +TtV(M) (3.1) 

where Am is the usual U(N) invariant Laplacian on the homogeneous space of hermitian 
matrices and the potential V(M) can actually explicitly depend on time t. 

The Schroedinger equation can be written in the form of a minimization principle: 

miny j d N2 MTr(^\d M ^(M)\ 2 + V(M)\q(M)\ 2 ) (3.2) 

To reduce this problem to the eigenvalues we use the U(N) symmetry of our model and look 
for a wave function \I/(M) transforming according to a certain irreducible representation 
RoiU(N): 

^(0+Mfi) = ^fi^(M) 
J 

where Or is a group element Q in representation R and /, J are the indices of the repre- 
sentation. Such a function may be decomposed as 

= (3-3) 
J 

Here if) R {xi, • • • , a; at) is a vector in the representation R. 

Near the unity element on the group space O ~ I + u we have Or ~ Pr + ■ oJijT^ 
where Pr is a projector (unity element) in the R space, ui is a small deviation from it and 
are the u(N) algebra generators. This gives: 



d . d A 1 d 



dMij dx k ^ x k - X m duOmj 

J m=l J 

and we finally obtain from ( |3.2j ) the following variational principle: 

J 7„ A J A -L A ™ 
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where all the quantities and operators with the subscript R are subjected to the corre- 
sponding matrix operations in the matrix space of representation. 
The Schroedinger equation now reads: 

- E A ~ 2 (*) ^^flST^*) - H T S T ^ X ) = ( E ~J2 v M)Mx) (3.5) 

k ' " i^j k 

It is useful to introduce a new function 4>r(x) — -^^''Pr(x) obeying the equation 

~ E j - E (^-^a ^) = ( E - E ^))^(*0 (3-6) 

Note that any translation u;^ — > u>ij + Sije does not change the wave function ^ R. That 
means that we are looking only for the states on which the condition is imposed 

T^ R = 0, k = l,---,N (3.7) 

On the first sight, we fulfilled our main task for the matrix quantum mechanics: we 
reduced it to an eigenvalue problem and are now dealing with only iV variables. But 
the Schroedinger equation ( |3.5| ) contains the Hamiltonian which is a matrix in the rep- 
resentation space acting on the wave function which is a vector in this space. For small 
representations whose Young tableaux contain << iV 2 boxes the problem is still solv- 
able in the large N limit (as we will demonstrate below). For a very interesting case of 
big representations (~ iV 2 boxes in the Young tableaux) the problem remains a serious 
challenge. 

In the simplest case of singlet representation (solved long ago in 0) the wave function 
is a scalar and the last term in the r.h.s. of the Schroedinger equation ( |3.6|) drops out. The 
problem appears to be equivalent to the quantum mechanical system of iV non-interacting 
fermions (due to the antisymmetry of (/>(x)) in a potential V(x). It was used in many 
applications, including the solution of the non-critical string theory in 1+1 dimensions 
0- 
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The next smallest representation is adjoint. The adjoint wave function satisfying the 
relation (|3.3|) should be a function of the type 



N-X 



a=0 

where the coefficients Ca possibly depend on the invariants (eigenvalues). If we denote 
4>adj( x i] x ) = 4>i{ x ) (depending of course on all iV £j) we can write the Schroedinger 
equation for the adjoint wave function in the form flTf] : 

+ V Mx) - ± £ = (3.8) 

One can see that the last term in the l.h.s. of this equation is ~ iV 2 smaller than the other 
terms and can be regarded as a small perturbation on the background of the free fermion 
solution of the singlet sector. 

For one of physically most interesting applications, the 1+1 dimensional string theory, 
we need to solve the model in the inverted oscillatorial potential V(M) = —M 2 . The model 
is unstable and one needs to specify the boundary conditions for big M's. Usually one 
considers the boundary conditions when the absolute value of any of the eigenvalues of 
M cannot exceed some maximum value A (a cut-off wall). In the case of the large iV 
limit one takes A ~ N and it happens that the spectrum density of the model depends 
in a very universal (logarithmic) way on A. In the singlet state the spectrum is that of 
N independent fermions (eigenvalues) in the same potential and the eigenfunctions are 
the Slater determinants of the parabolic cylinder functions (see the review in p| . In the 
non-singlet sectors the eigenvalues start interacting and obey a more complicated statistics 



corresponding to the symmetry of the Young tableau of representation (see the review [|19 



for the details). Although the problem is clearly integrable the spectrum of the non-singlet 
sectors of the inverted matrix quantum oscillator is still unknown (for the large N estimates 



of the mass gap of adjoint representation see Hl7|,p0| and pl| ). 

It was conjectured in [^Of and shown in [21| that the adjoint representation describe the 
vortex anti- vortex sector in the 1+1 dimensional string theory with one compact dimension. 
Higher representation describe higher numbers of vortex anti-vortex pairs (corresponding 
to the number of boxes in the Young tableau of the representation). 



4. Character expansion and new solvable (multi) matrix models 

The group character expansion has shown its power in the lattice gauge theory long 
time ago, starting from the work of A. Migdal [p2 |. 



The character expansion method proposed in the papers [^3[- |24| and inspired by the 



result of paper |25[ is the most general approach for the reduction of the number of degrees 
of freedom from ~ N 2 to ~ N in a new big class of (multi) matrix integrals. The matrix 
integral considered in these papers looks as follows: 

N 2 i\ r — .„ r itl, ji /f2 



Z = J cT Mexp[-TrM 2 + TrV {AM)] (4.1) 

where V(y) = '^2k > 2^ky k is an arbitrary potential and A is an arbitrary hermitian matrix 
(which can be taken diagonal without a loss of generality). We again diagonalize the matrix 

M as 

m = n + xn (4.2) 

The integral over the U(N) variable O looks difficult to do directly since the Itzykson-Zuber 
formula ( |2.7|) seems to be of little use here. Instead of it let us expand exp[TW(AM)] as 
an invariant function of the variable AM in terms of the characters xr{AM) of irreducible 
representations R of the GL(N) group: 

exp[TW(AM)] = fRXR(AM) (4.3) 

R 

where the coefficients fn are the functions of iV highest weight components of a represen- 
tation 

R = {0 < rriN < rriN-i < • - ■ < rti\ < oo}. 

The sum ^2 R is nothing but the sum over N ordered integers. They can be calculated due 
to the orthogonality of characters as the following unitary integrals: 

f R = J [dn} u(N) exp[TW(fi)] X K(0+) (4.4) 



Note that this integral can be represented as an explicit integrals only over the Cartan 
subgroup O = {e 1 ^ 1 , • • • ,e tuJN } and thus contains only N integration variables. We have 
TrV(0) = J2 k V(e^) and [dn] u{N) -> Yl k dd k Yl^j sin 2 ^=2i. Now if we plug ([Op into 
( [OP we realize that the decomposition (|4.2|) is actually useful and we can integrate over 
Q using the following orthogonality relation between matrix elements of representation R: 

f [dn] u{N) x R (An+xn) = -J— Xr (a) Xr (x) (4.5) 

where dirriR is the dimension of a representation R. We see that we achieved our main 
goal: due to the formulas ([0|), (|4.4j ) and ( f4.5| ) we reduced the original matrix integral 
( |4.1| ) to an integral over only N eigenvalues xi, ■ ■ ■ , xn of the matrix M and the sum over 
N highest weight components mi, • • • , toat- In the large N limit, if we scale appropriately 
the constants in the potential V(M), the sums over m's can be replaced by integrals and 
we can again apply the saddle point approximation in all 2N integration variables. To get 
explicitly the right large iV scaling of the couplings one usually changes e v — > e NV . Then 
the effective action at the saddle point is always of the order 1/N 2 and the new couplings 
of the potential V can be kept finite in this limit. 



As was shown in [£5| (see also P3|), the integral over x\, ■ ■ ■ , xn can be calculated 
exactly and the remaining sum over strictly ordered nonnegative integers hi = —rrii +N — i 
(shifted highest weights) reads: 

V\(h e - l)\\h°U 

Z = E "n (fe « Jw" x*Wx*lt) (4.6) 

h 1 <h 2 <---<h N 11V ^ 

where {h e } and {h°} are the collections of even and odd integers hk (their number is 
equal). Only the representations with equal amounts of even and odd /t'es contribute to 
( fO| ) . The products in the numerator go over all even and odd /i'es and the product in the 
denumenator goes over all couples h e ,h° . xn(t) is a character of the coupling constants tu 
written in the Schur form: 

X R = det ij P hi - j (t) (4.7) 
9 



, . ^ * „, fc fc . 

So in the large iV limit we have to do the saddle point calculation only with respect 
to iV summation variables hi, ■ ■ ■ , hjy. 

The details of these formulas can be found in [2S|, R3]-[23]. One can also find in 



these papers the geometrical interpretation of the integral (|4.1|) in terms of the so called 
dually weighted planar graphs. It gives the generating function of planar graphs where 
both vertices and faces are weighted by the generating parameters depending on their 



orders. In ||23|1 - ||24|| one can find the solutions of some combinatorial problems related to 
the enumeration of planar graphs which were possible only due to the power of the character 
expansion method. The particular solutions of the saddle point equations could be very 
tricky but it is already a "classical" problem of solution of various integral equations rather 
than a "quantum" problem of functional integration over infinite matrices. In that sense 
this model is solvable. 

It is obvious that there exist many ways to generalize the model (|4.1|) to other matrix 
integrals. An immediate generalization is to substitute the TrM 2 term in Q4.1| ) by an 
arbitrary function W(M). In that case we cannot calculate explicitly the coefficients fn 
(except when W is a monomial: W(M) = M k ) but we still get an explicit integral over 
3N variables Xi, u>i and mj. So the model is again solvable. 

Another solvable matrix model of this kind involving general complex matrices was 



proposed and investigated in [26] , [g7f . Its free energy gives a generating functional count- 



ing branched coverings of two dimensional surfaces. 

The most general solvable two matrix model reads as 

Z = J d N2 Ad N2 BexpNTr(U(AB) + V(A)+W(B)) (4.8) 

where U,V and W are arbitrary functions. The way to reduce it to ~ N degrees of freedom 
is again to expand in characters 

exp[TW( AB)] = urXr(AB) (4.9) 

R 
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diagonalize the matrices A and B and integrate over the U(N) variable between them by 
means of (O) . In a particular case 



Z= f d N2 Ad N2 BexpNTr(-(A 2 +B 2 )-j(A 4 + B 4 )-P(AB) 2 ) (4.10) 
J 2 4 2 

the model describes a special trajectory of the 8-vertex model on random graphs. It was 
completely solved in E3] . Again it was possible, using character orthogonality relations, to 
integrate over the relative angle between A and B; this leads to separation into one- matrix 
integrals: 

Z(a,fi ~ ^(iV/3/2)^/ 2 c {M [P { , } (a)] 2 (4.11) 
{h} 

where c^} is a coefficient: 

1 

C{h} ~ n^mn^hr* - hf d ) 



and S^jfa) is a one-matrix integral 

%}(«) = J d N M X{h} (M)expN 



--trM 2 + -trM 4 
2 4 



(4.12) 



which appears squared in ( |4.11| ) because the contributions from the two matrices A and 
B are identical. 

Now we can reduce the calculation of the one-matrix integral S{h} to eigenvalue inte- 
grations: 



(4.13) 



where A(A) = det (^" J j = n i<& (A,- - X k ). 

Now we are left only with degrees of freedom and the action of the order A 2 , so 
the integration is reduced to the saddle point calculation with respect to the eigenvalues 



Afc (see P|| for the details). 
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We can immediately propose some solvable generalization of the general two matrix 
model (|4.8|) to a multi-matrix chain: 

f Q Q Q-1 

Z= l[d N2 M q expTv(j2v q (M q ) + J2W(M P _ 1 M P )) (4.14) 

J 9 =1 g=l p=l 

where V and W are arbitrary functions. 

Another interesting model solvable by the character expansion method can be written 
in the following general form: 

, Q Q 
Z = / II rf7V Hexp (TrV q (M q ) + V( J] M p )) (4.15) 

J g=l p=l 

It can be solved by character expansion with respect to the last factor by means of the 
formulas ( |4.3p ,( JOD and the multiple application of the formula (|4.5[ ) by induction. 



The results in terms of the sum over N highest weight components of the representa- 
tions R reads: 



where 



S W = j dN2M X{h}(M)expTrV q (M) (4.17) 



The last integral can be immediately reduced to the integrations of a type ( |4.13| ) over 
eigenvalues of the matrix M. 



5. Comments and unsolved problems 

A few comments are in order: 

1. The non-singlet sectors in the matrix quantum mechanics ( |3.1j ) can be effectively 
studied for the oscillatorial potential V(M) = M 2 . In this case the Hamilton-Ian is a 
collection of A^ 2 independent oscillators represented by the matrix elements of M. The 
spectrum of Hamilton-Ian of this model in a given irreducible representation of U(N), is 
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encoded into the partition functions Zf>(q) for finite inverse temperature (3 (where q = 
in a given representation R. The effective way to study Zn(q) can be found in papers [^1 
or 



29]. 



2. The character expansion is nothing but the Fourier expansion on a group manifold. 
As trivial as it looks for us now the Fourier transform was always a powerful method of 
solving problems using their symmetries. Many of the matrix models presented in the 
previous section and solved by this method seemed hopeless just a few years ago. 

3. One of the interesting and not well studied questions is how to classify all the 
matrix integrals which can be reduced from ~ N 2 to ~ N integrals or sums by use of the 
character expansion. 

4. In many physically interesting cases we don't need a general form of potentials 
V(M) or W(M) mentioned through this paper. For example, as it was mentioned, to 
study the universal behavior of the large N matrix quantum mechanics near the instability 
point we need to know only the solution in the vicinity of a quadratic top of the potential 
V(M) ~ — Tr(M - M ) 2 . The rest of the potential V(M) has small influence on the 
behavior of the eigenvalues and serves only as a U(N) invariant cutoff wall. It simplifies 
greatly the problem. For instance, all the applications in string theory, two-dimensional 
quantum gravity and most of statistical-mechanical applications need only the analyses of 
the vicinity of such critical points. The lesson to draw from it is that for some physically 
most interesting regimes the seemingly hopeless matrix integrals become not so hopeless 
and look "almost Gaussian". May be a general method of the investigation of these 
instability points can be worked out. 

5. Another question is related to the integrability properties of sums and integrals after 
such reduction. The partition functions of some of them (such as the old one matrix and 
two matrix models) are known to be r-functions of some integrable hierarchies of classical 
differential or difference equations, like Toda hierarchy ( |4.15| ), ||29|| or KP hierarchy |30| (see 
31|]for a good introduction). But many others, like the model (^4.1| ), cannot be represented 
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by free fermions. On the other hand, the Itzykson-DiFrancesco formula ( f4.6| ) suggests 
that it might exist some interacting fermion representation of the partition function of the 



model of dually weighted graphs (|4.1 



The method of character expansion as well as all other methods of calculation of the 
large iV matrix integrals presented here represent just another refining and generalization 
of the usual method of reduction to the matrix eigenvalues invented long ago by Dyson. 
Its range of applicability is quite limited although it includes quite a few important matrix 
integrals known from physics and mathematics. Many more interesting matrix integrals 
look not hopeless for the investigation in the large iV limit. The search for new tricks 
of integration over matrices is a fascinating and potentially extremely rewarding research 
direction. 
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